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ABSTRACT 

We report FUSE far-UV spectroscopy of the prototypical dwarf starburst galaxy 
NGC 1705. These data allow us for the first time to directly probe the coronal-phase 
(T = few xlO 5 K) gas that may dominate the radiative cooling of the supernova-heated 
ISM and thereby determine the dynamical evolution of the starburst-driven outflows in 
dwarf galaxies. We detect a broad (~100 km s _1 FWHM) and blueshifted (Av = 77 km 
s _1 ) Oy/A1032 absorption-line arising in the previously-known galactic outflow. The 
mass and kinetic energy in the outflow we detect is dominated by the warm (T ~ 10 4 
K) photoionized gas which is also seen through its optical line-emission. The kinematics 
of this warm gas are compatible with a simple model of the adiabatic expansion of a 
superbubble driven by the collective effect of the kinetic energy supplied by supernovae 
in the starburst. However, the observed properties of the OVI absorption in NGC 1705 
are not consistent with the simple superbubble model, in which the OVI would arise 
in a conductive interface inside the superbubble's outer shell. The relative outflow 
speed of the OVI is too high and the observed column density (logNovi = 14.3) is 
much too large. We argue that the superbubble has begun to blow out of the ISM of 
NGC 1705. During this blow-out phase the superbubble shell accelerates and fragments. 
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The resulting hydro dynamical interaction as hot outrushing gas flows between the cool 
shell fragments will create intermediate-temperature coronal gas that can produce the 
observed OVI absorption. For the observed flow speed of ~ 10 2 km s _1 , the observed 
OVI column density is just what is expected for gas that has been heated and which 
then cools radiatively. Assuming that the coronal-phase gas is in rough pressure balance 
with the warm photoionized gas, we estimate a cooling rate of-order ~0.1 M & per year 
and ~10 39 erg s _1 in the coronal gas. The latter represents <10% of the supernova 
heating rate. Independent of the assumed pressure, the lack of observed redshifted 
OVI emission from the backside of the outflow leads to upper limits on the cooling 
rate of < 20% of the supernova heating rate. Since the X-ray luminosity of NGC 1705 
is negligible, we conclude that radiative-losses are insignificant in the outflow. The 
outflow should therefore be able to fully blow-out of the ISM of NGC 1705 and vent its 
metals and kinetic energy. This process has potentially important implications for the 
evolution of dwarf galaxies and the IGM. 

Subject headings: galaxies: individual: NGC 1705 - galaxies: starburst - galaxies: dwarf 
- galaxies: kinematics and dynamics - galaxies: halos - galaxies: ISM 

1. Introduction 

Local starburst galaxies are excellent local laboratories to study the physics of galaxy building 
(cf. Heckman 1998). It is now clear that the deposition of mechanical and thermal energy by 
multiple supernovae in starbursts leads to a global outflow of metal-enriched gas. These outflows 
are called "superbubbles" during their early dynamical evolution, and "superwinds" after they 
blow out of the galaxy's ISM (e.g., Heckman 2000 and references therein). Such flows are expected 
to play an especially important role in the evolution of dwarf galaxies, whose relatively shallow 
potential wells make them susceptible to wind-driven loss of gas and newly-created metals (e.g., 
Dekel & Silk 1986; Martin 1999). Since similar outflows appear to be common in high-redshift 
galaxies as well (Pettini et al. 1998, 2000; Tenorio-Tagle et al. 1999) they are perhaps the most 
plausible mechanism by which the mass-metallicity relation in galactic spheroids was established 
(e.g., Lynden-Bell 1992) and the intergalactic medium was heated and chemically-enriched (e.g., 
Ponman, Cannon, & Navarro 1999; Gibson, Loewenstein, & Mushotzky 1997). 

One of the major uncertainties concerning starburst-driven outflows is the importance of ra- 
diative cooling: what fraction of the kinetic energy supplied by supernovae is carried out in the 
flow rather than being radiated away? While the available X-ray data and models establish that 
radiative losses from hot gas (T > 10 6 K) are not severe (see Strickland & Stevens 2000), up until 
now there has been no direct observational probe of the coronal-phase gas (T = 10 5 to 10 6 K) that 
could dominate the radiative cooling. 

The recent launch of the Far Ultraviolet Spectroscopic Explorer (FUSE - Moos et al. 
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2000) provides access to the best probe of rapidly-cooling coronal gas in starburst outflows: the 
OVIAA1032,1038 doublet. Accordingly, we have obtained FUSE spectra of a small sample of the 
nearest and brightest starbursts. In this paper, we present the first FUSE detection of OVI in 
the galaxy NGC 1705. 

NGC 1705 is an ideal test-case for the study of coronal-phase gas. This nearby (D = 6.2 Mpc) 
dwarf starburst galaxy has the second highest vacuum-ultraviolet flux of any starburst galaxy in the 
extensive compilation of IUE spectra in Kinney et al. (1993). The detailed investigation by Meurer 
et al. (1992) established NGC 1705 as a prototypical example of a dwarf starburst undergoing 
mass-loss. They were able to delineate a kpc-scale fragmented ellipsoidal shell of emission-line 
gas that was expanding at roughly 50 km s _1 along our line-of-sight. They also showed that the 
population of supernovae in the young super star cluster (NGC 1705-1) was energetically-sufficient 
to drive this flow. The expulsive nature of the flow was later confirmed by HST observations that 
showed that the ultraviolet interstellar absorption-lines towards NGC 1705-1 were blueshifted by 
70 to 80 km s _1 relative to the galaxy systemic velocity (Heckman & Leitherer 1997; Sahu & Blades 
1997; Sahu 1998). Hensler et al. (1998) reported the detection of soft X-ray emission, presumably 
from hot gas inside the expanding emission-line nebula. 



2. Observations & Data Reduction 

Two FUSE observations of NGC 1705 (a 20 oo = 04 /l 54 m 13.48 s , 5 20 oo = -53°2l m 39.4 s ; l n = 
261.0788°, bn = —38.7428°) were obtained on 4-5 February 2000. The central super star cluster 
NGC 1705-1 was centered in the large (LWRS, 30" x 30") aperture of the LiFl (guiding) channel for 
each observation by the standard guide-star acquisition procedure. The two observations resulted 
in 13 exposures totaling 21.3 ksec of on-target exposure time. Approximately 95% of the observing 
time occurred during orbital night, which greatly reduced the amount of terrestrial O I and N I 
airglow entering the FUSE apertures. Flux was recorded through the LWRS apertures in both 
long wavelength (LiF, ~ 1000 — 1187 A) channels and both short wavelength (SiC, ~ 900 — 1100 A) 
channels. The astigmatic heights of the LiFl spectra near 1030 A were roughly 1/3 of the aperture 
width, consistent with the compact photometric structure of NGC 1705. 4 The data are preserved 
in the FUSE archive with observation identifications A0460102 and A0460103. 

The raw time-tagged photon event lists for each exposure were processed with the standard 
FUSE calibration software (CALFUSE vl.7.5) available at the Johns Hopkins University as of 



4 HST images at 2200A show that roughly 40% of the light that would be admitted into the FUSE LWRS aperture 
would come from the single brightest star cluster NGC 1705-1, and 70% of the light would come from the region 
within a radius of 2.8 arcsec of NGC 1705-1 (Meurer et al 1995). A comparison of the HST GHRS (1.7 arcsec 
aperture) and IUE (10 by 20 arcsec aperture) spectra implies that these fractions should not be strongly wavelength 
dependent. Thus, while the FUSE data sample the entire inner region of NGC 1705 (projected size of roughly 900 
by 900 pc), the sightlines within about 80 pc of NGC 1705-1 are very strongly weighted. 
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August 2000. The lists were screened for valid data, and corrections for geometric distortions, 
spectral motions, and Doppler shifts were applied (see Sahnow et al. 2000). The data were not 
affected by the detector event bursts that plagued many of the early FUSE spectra. The 13 
individual calibrated extracted spectra for each channel were cross-correlated, shifted to remove 
residual velocity offsets due to image motion in the apertures, and combined to a produce composite 
spectrum for each channel. These composite spectra in the 1000-1070 A region were then compared, 
and any remaining velocity offsets were removed by referencing to the LiFl channel data. The 
wavelengths were put into the Local Standard of Rest (LSR) reference frame by requiring that 
the Galactic ISM lines fall at ~ +20 km s _1 , the approximate velocity of the Milky Way lines 
as observed at longer wavelengths with the Space Telescope Imaging Spectrograph (STIS) on the 
HubbleSpaceTelescope. These new STIS data will be presented in future papers (Sembach et al. 
2001; Heckman et al 2001). The FUSE data have a wavelength uncertainty of ~ 6 km s _1 (la). 
The spectrum of NGC 1705 is shown in Figure 1. 

In this article, we do not combine the data from all four channels because the instrumental 
resolution and sensitivity changes as a function of wavelength. Therefore, when measuring the 
strengths of absorption features, we compare the individual measurements (W^) for the two chan- 
nels having the highest sensitivity at the wavelengths of the absorption features of interest (usually 
LiFl and LiF2 for A > 1000 A, or SiCl and SiC2 for A < 1000 A). The integrated equivalent widths 
derived from separate channels generally agree very well. Table 1 contains these measurements for 
selected Milky Way and NGC 1705 absorption features observed along the sight line. The data 
have a velocity resolution of ~ 30 km s _1 and S/N ratios of 16, 13, 10, and 9 per resolution element 
at 1032 A in the LiFl, LiF2, SiCl, and SiC2 channels, respectively. 

3. Results 

We confirm the principal conclusion of Sahu &: Blades (1997) and Sahu (1998): there are three 
separate systems of interstellar absorption-lines along the sight-line to NGC 1705 (Table 1; Figure 
2). These arise in the Milky Way, the periphery of the high velocity cloud HVC 487, and the 
interstellar medium of NGC 1705. In the present paper, we will focus our analysis on NGC 1705, 
but briefly summarize our results on HVC 487. 

3.1. HVC 487 

As pointed out by Sahu & Blades (1997), HVC 487 is located 2° away from the NGC 1705 
sight-line. Sahu (1998) proposes that HVC 487 is associated with the Magellanic Stream, so the 
implied impact parameter would be roughly 2 kpc. The properties of the absorption associated 
with HVC 487 will be discussed in detail in a future paper which will combine the FUSE data with 
echelle spectra taken with STIS on HST (Sembach et al. 2001). Here, we note only that HVC 487 
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Fig. 1. — Overview of our FUSE spectrum on NGC 1705. These data have been binned into 0.04 A 
samples for clarity. Data from sing le segments are shown (SiC2A: 910-1005 A, LiFIA: 995-1085 A, 
SiC2B: 1075-1090 A, LiF2A: 1090-1182 A). The convergence of the Lyman series of hydrogen is 
evident, as are numerous strong interstellar lines arising in both the Milky Way and NGC 1705. 
The upward spike near 1026 A is due to Lyf3 airglow. 
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Fig. 2. — FUSE data for the four detector segments covering the 1018-1045 A spectral region. 

These data have been binned into 0.04 A samples. Prominent absorption lines of atomic species 

in the ISM of the Milky Way, HVC 487, and NGC 1705 are marked above the LiFIA spectrum. 

The HVC is seen only in the lines of Ly/3, OVI, and CII in the spectral region shown. HI Ly(3 

airglow present near zero velocity can be used to judge the spectral resolution for a source filling 

the LWRS apertures (FWHM ~ 120 km s^ 1 ). The data have a (point source) spectral resolution 

of ~ 30 km s" 1 and S/N ~ 20-25 per resolution element in the wavelength region shown. The 

dashed line over-plotted on the top spectrum indicates the continuum adopted for our analyses of 

the absorption lines. 
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is strongly detected in absorption in the OVIX1032, CIIIX977, CIIX1036, and Lyman-series lines. 
Assuming that the OVIXW32 line is optically-thin, the implied OVI column is 2.0±0.3 x 10 14 
cm" 2 . This is typical of other HVC sight-lines studied with FUSE (Sembach et al. 2000). The 
Lyman series, CIIIA977, and CIIA1036 absorption-lines from HVC 487 are at vlsr = 270±15 km 
s -1 . These are close to the HIX2lcm velocity of HVC 487 of 232±14 km s" 1 (Sahu 1998). In 
contrast, the OV/A1032 line is centered at vlsr = 326±10 km s _1 , an offset of 94±17 km s~ x from 
the HIX21cm velocity. The OVI line is very broad, with an observed FWHM = 100±15 km s~ x . 
This breadth is comparable to the sound speed in coronal gas where OVI would be abundant, and 
is similar to those of other OVI HVCs observed with FUSE (Sembach et al. 2000). 

3.2. NGC 1705 

We expect our FUSE data to probe four phases of interstellar gas in NGC 1705: molecular 
gas, neutral atomic gas, warm gas photoionized by hot stars, and collisionally-heated coronal gas. 
The neutral atomic gas is traced by species with ionization potentials of creation x < 1 Ryd, the 
warm gas by ions with % = 1 to 4 Ryd (up to the Hell edge), and the coronal gas by ions with 
x >4 Ryd. 

No molecular hydrogen was detected in any of the first three J levels (0-2), and the limit on 
the total H2 column density is logNn2 < 14.6 (corresponding to a molecular hydrogen fraction of 
J'h2 = 2Nh2/(Nhi + 2Nh2) < 6 x 10~ 6 - based on the observed HI column given below). Similar 
results have been reported for the metal-poor dwarf starburst IZwl8 (Vidal-Madjar et al. 2000). 
The intrinsic reddening in NGC 1705 is very small: Heckman et al. (1998) give E{B — V) « 
0.1. In the particular case of NGC 1705, the low values for Nh2 and ju2 are quite consistent with 
Galactic sightlines with similarly- low E(B — V) and Nhi+H2 (Savage et al. 1977). 

Lines tracing the other three phases are well-detected (Table 1; Figures 1 and 3). In particular, 
the Oy/A1032 line is independently detected at the IOct level in both the LiFl and LiF2 detectors 
(Figure 2). The weaker member of the doublet (OF/A1038) is blended with the much stronger 
Galactic OIA1039 line and is not convincingly detected. 

3.2.1. Kinematics 

We begin by considering the kinematic properties of the gas. In what follows, we will adopt 
an LSR systemic velocity for NGC 1705 of v sys = 622 km s _1 , which is derived from the HIA21cm 
rotation curve (Meurer et al. 1998). 

The mean radial velocity of the neutral gas lines listed in Table 1 is 590±11 km s _1 as deter- 
mined by fitting Gaussian components to the observed absorption features. This implies a blue-shift 
of 32 km s _1 relative to v sys , which is significantly larger than the uncertainty 11 km s" 1 derived 
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Fig. 3. — Continuum-normalized absorption line profiles versus LSR velocity for selected absorption 
lines in the FUSE bandpass. The profiles shown are from the LiFl channel, except for CIII A977.0, 
OIA988.8, and NIIA1084.0, which are from the SiC2 channel. Data of roughly comparable quality 
exists for the LiF2 and SiCl channels (see Table 1). The identifications of the lines are listed 
underneath each spectrum. Additional lines arising in the ISM or within NGC 1705 within the 
velocity range shown are also indicated (without wavelengths). In some cases, lines from various 
species blend and cannot be resolved. The vertical dashed lines indicate the LSR velocities of the 
ISM (20 km s" 1 ), HVC 487 (~ 305 km s" 1 compared to the i27A21cm emission velocity of 232 
km s" 1 ), and the systemic velocity of NGC 1705 (622 km s" 1 ). Note that the NGC 1705 absorption 
lines are substantially blue-shifted relative to the systemic velocity of the galaxy. 
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from the standard deviation of the line centroids. In contrast, the coronal gas traced by OVI 
is significantly more blueshifted: vovi = 545+10 km s _1 , or v — v sys = -77 km s _1 . The warm 
photoionized gas appears to have intermediate velocities with the mean line centroid at 569+10 km 
s _1 (v — v sys = -53 km s -1 ) 5 . 

Further evidence for differing dynamics between the neutral and coronal gas comes from the 
line widths (see Figure 3) . The weak (unsaturated) lines from the neutral phase have observed lines 
widths (FWHM) of roughly 70 km s -1 , while the Oy/A1032 line has a broader structure (FWHM 
~ 100+10 km s _1 ). The lines from the photoionized gas are roughly the same width as the OVI 
line. 

The strongly saturated CIIIX977 and CI/A1036 lines allow us to probe the low-column density 
gas at extreme radial velocities. Relative to v sys , the centroid of the CIII line is blueshifted by 
95+20 km s _1 . The blueward edge of the line is blended with the HVC feature, but absorption 
extends to at least -260 km s -1 with respect to v sys . Relative to v sys , the CT/A1036 line centroid is 
blueshifted by 40+20 km s _1 km s _1 and absorption is present from at least -180 km s -1 (where the 
line is blended with the HVC feature) to roughly +70 km s _1 . The HIA21cm maps of NGC 1705 
published by Meurer et al (1998) show emission from -94 to +86 km s ^ relative to v S y S . The 
velocity range of the redshifted gas seen in C//A1036 is thus consistent with the kinematics of the 
normal neutral (turbulent?) ISM, but an additional component of highly blueshifted (outflowing) 
neutral gas is also present. 



3.2.2. Column Densities and Abundances 

We have determined column densities by converting the observed absorption profiles into op- 
tical depth profiles and integrating over velocity (see Savage & Sembach 1991). We measured only 
the weaker (and thus, more-optically-thin) metal lines with equivalent widths W < 250 mA. The 
results are given in Table 2. 

For the lines in the neutral phase, we estimate the following column densities (logiV): NI (14.0), 
OI (15.6), Sill (14.6), Arl (13.5), and Fell (14.5). Heckman k Leitherer (1997) estimated that 
Nhi = 1.5xl0 20 cm -2 column towards NGC 1705-1 based on fitting the red side of the damped 
Lya profile in the HST GHRS data. We have fit the Ly/3 and higher-order Lyman series lines in 
our FUSE spectrum, and find Njji = 2+1 xlO 20 cm~ 2 (6 = 35+5 km s _1 ). We adopt the mean 
of these two estimates: logNm = 20.2+0.2. Since these species above should be the dominant 
state of their respective element in the neutral gas, we can estimate the following abundances 
(Anders & Grevesse 1989): [N/H]=-2.2, [0/H]=-1.5, [Si/H]=-l.l, [Ar/H]=-1.3, and [Fe/H]=-1.2. 
The absolute abundances are uncertain by 0.3 dex. With the exception of nitrogen, the abundances 
are marginally consistent with the metallicity of [O/H]=-0.9 derived for the nebular emission-line 



5 We exclude the highly saturated and very broad CIIIX977 line, which is discussed separately 
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gas (Heckman et al. 1998). This is in agreement with expectations that we see ambient interstellar 
gas, rather than the much hotter gas that is polluted by recent supernova ejecta. 

We confirm Sahu's (1998) finding of a low relative N abundance, based on the NI column 
density. As Sahu points out, the Nil ion can also make a significant contribution to the total N 
column in the HI phase (see Sofia &: Jenkins 1998). We estimate logNwu = 14.7. Unless most of 
this is associated with the HI phase, N remains selectively under abundant. In fact, most of the 
observed Nil must be associated with the warm photoionized gas that dominates the total gas 
column along our line-of-sight (see below). This is implied by the kinematics of the Nil A1084 line, 
which has a significantly larger blueshift and FWHM than the lines arising in the neutral phase, 
but agrees with the other lines arisiing in the warm photoionized gas (Table 1; and see above). A 
low nitrogen abundance is typical of low-metallicity gas in dwarf galaxies, and is a consequence of 
a partly secondary nucleosynthetic origin for N (e.g. van Zee, Haynes, & Salzer 1997). We do not 
confirm the low relative Fe abundance noted by Sahu &: Blades (1997). Finally, we see no evidence 
for the systematic gas-phase depletion of refractory elements (e.g. compare Fe to Ar). The neutral 
gas in the outflow thus appears largely dust-free (e.g., Savage &: Sembach 1996), consistent with the 
lack of detectable reddening in the UV spectrum of NGC 1705 (e.g., Meurer, Heckman, & Calzetti 
1999). 

The total column density in the warm photoionized gas is more uncertain. The weak S7iTA1012 
and S , /UA1063 lines imply column densities logNsni = 15.0 and logNsiv = 14.4. The SII ion 
can also be adundant in photoionized gas. While the HST GHRS spectrum covers a somewhat 
different sight-line (the 1.7x1.7 arcsec Large Science Aperture was centered on NGC 1705-1), the 
SII column density logNsn = 15.0 from Sahu & Blades (1997) is a useful additional constraint. 
The total implied column of ionized S is logAs = 15.3. For an assumed value [S/H]=-0.9 (based 
on nebular emission-line metallicity) , the total HII column in the photoionized gas is logNnii = 
20.9, or almost an order-of-magnitude greater than the HI column. 

The observed OVIXW32 line yields an OVI column density logNovi = 14.3. Deriving the 
total column density of the coronal phase gas requires uncertain assumptions. For gas in collisional 
ionization equilibrium, OVI reaches its peak relative abundance (~20% of the total O abundance) 
at T ~3xl0 5 K (Sutherland & Dopita 1993). If we assume that [O/H]=-0.9 in the coronal phase 
gas (similar to the value in the optical emission-line gas), the implied minimum total HII column 
in the coronal gas is then logiV cor = 19.0. While uncertain, we note that this is a substantially 
smaller column than seen in the cooler gas. 

3.3. Limits on OVI Emission 

In addition to the blueshifted OVI absorption-line produced by coronal gas on the frontside 
of the outflow, redshifted OVI emission from the backside of the outflow must also be present. 
We do not detect such emission, but can set an interesting upper limit to its intensity. Combining 
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our LiFl and LiF2 spectra, we obtain a 3cr upper limit on the flux of the redshifted OV/A1032 
emission-line of < 7.5xl0" 15 erg cm" 2 s _1 . This limit assumes that the line would have a breadth 
similar to the blueshifted absorption-line (~ 100 km s" 1 ). 

We first correct this flux for forground dust extinction. The Galactic HI column towards 
NGC 1705 is 1.3xl0 20 cm" 2 , based on both the damped Lya profile (Sahu 1998) and radio A21cm 
observations (reported in Hensler et al. 1998). Adopting the standard extinction curve of Mathis 
(1990), the implied extinction at 1032A is 0.3 magnitudes. The data discussed in section 3.2.2 
imply that there is a negligible amount of dust extinction intrinsic to NGC 1705. 

In the discussion below, we will be interested in the total OVI luminosity of NGC 1705. To 
derive a global upper limit, we make a simple aperture correction to our observed limit. The portion 
of the backside of the NGC 1705 outflow included in our FUSE aperture represents about 15% of 
the total surface area of the outflow (front and back sides), as mapped by Meurer et al. (1992). 
Thus, the global extinction-corrected OV/A1032 flux from NGC 1705 is < 6.5 xlO -14 erg cm -2 
s" 1 , and the corresponding luminosity is < 3xl0 38 erg s" 1 . 

4. Discussion 

4.1. A Simple Model 

The simplest physical model that can be compared to our data is an adiabatic expanding 
superbubble whose expansion is driven by the energy supplied by multiple supernovae (Weaver et 
al. 1977; Koo & McKee 1992). This is a natural model to use, since it appears to be a good 
first-order description of the properties of the Ha emission line nebula in NGC 1705 (Meurer et al. 
1992; Marlowe et al. 1995). 

In this simple model, there are six concentric zones. From inside-out these are: 1) An innermost 
region inside which energy is injected by supernovae (the starburst). 2) A region of supersonic flow 
fed by the hot gas created in zone 1. 3) A region of hot gas (material from zone 2 that has passed 
through an internal shock). 4) A conductive interface of intermediate-temperature gas created as 
the hot gas in zone 3 heats the relatively cool, dense gas in zone 5. 5) A thin dense shell of ambient 
gas that has been swept-up (shocked and then radiatively cooled) as the "piston" of hot gas in zone 
3 expands into the ISM. Depending on the available flux of ionizing radiation from the starburst, 
this shell may be partially or fully photoionized. 6) The undisturbed ambient ISM. 

The dynamical evolution of a superbubble in the type of plane-parallel ISM appropriate to a 
starburst has been extensively discussed (see MacLow &: Ferrara 1999 and references therein). Once 
the radius of the superbubble is several times the vertical scale-height of the ISM, its expansion will 
accelerate and Rayleigh- Taylor instabilities will cause the outer shell (zone 5) to fragment. This 
will allow the hot gas from the interior to escape from the ruptured superbubble and flow out into 
the galactic halo. This "blow-out" or " break-out" stage marks the transition from superbubble to 
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superwind. 

4.2. Some Simple Inferences 

Before considering the superbubble model in more detail, it is worthwhile to make some rough 
estimates of the mass and kinetic energy in the gas and compare these to expectations. To convert 
the observed column densities and outflow velocities into masses and kinetic energy, we will assume 
the following idealized model for the NGC 1705 outflow. Meurer et al. (1992) show that the 
morphology and kinematics of the emission-line nebula can be described as the expansion of a hollow 
prolate ellipsoid with a semi-major axis of 1.5 kpc and semi-minor axes of 0.5 kpc. Multiplying the 
surface area of this ellipsoid by our estimated column densities (see above) implies total gas masses 
of 5x10 , lxlO 7 , and 6xl0 5 Mq respectively in the warm photoionized gas, the neutral gas, and 
the coronal gas. The outflow speeds then yield corresponding kinetic energies of 1.6xl0 54 , lxlO 53 , 
and 3xl0 52 ergs. 

The supernovae in the NGC 1705 starburst will supply kinetic energy at a mean rate of 
roughly 2xl0 40 erg s _1 (Marlowe, Meurer, Sz Heckman 1999). The dynamical age of the expanding 
emission-line nebula is 10 to 15 Myr (Marlowe et al. 1995; Meurer et al. 1992), so the total amount 
of kinetic energy supplied by the starburst during this time is ~ 1 x 10 55 ergs. In the standard 
adiabatic superbubble model, the kinetic energy of the swept-up shell (zone 5) will be about 19% 
of the total injected energy (e.g. MacLow & McCray 1988). Given the nature of our estimates, the 
rough agreement between the observed and available/predicted kinetic energies is gratifying. 

The amount of mass returned directly by supernovae and stellar winds during the last 10-15 
Myr will only be ~10 5 Mq (Leitherer & Heckman 1995). This is tiny compared to shell's mass (see 
above), consistent with the basic superbubble model in which the shell is swept-up ambient ISM. 
In this case, the total gas column density seen by FUSE (10 21 cm~ 2 ) can be no larger than the 
total column density of the ISM prior to the expansion of the superbubble. The HIA21cm maps in 
Meurer et al. (1998) show that the mean column density in the inner part of NGC 1705 is about 
2xl0 21 cm -2 , consistent with this requirement. 

The similarity between the total gas column of outflowing gas measured by FUSE and the 
typical HI column in NGC 1705 implies that the dimensions of the superbubble must be at least 
comparable to the characteristic thickness of the ISM. This suggests that the superbubble is early in 
the blow-out stage of it dynamical evolution. Such an inference is also consistent with the patchy, 
filamentary morphology of the superbubble's Ha emission (Meurer et al. 1992; Marlowe et al. 
1995), which suggests that the superbubble shell has begun to accelerate and fragment (see above). 

The absorbing HI column in NGC 1705 (logNni = 20.2±0.2) is roughly an order-of-magnitude 
smaller than the column observed in emission at A21cm in the same region (with the ~30 arcsec 
ACTA beam - Meurer et al. 1998). This implies that roughly 90% of the HI seen in the radio 
map lies behind the dominant source of far-UV light (the region around the super star cluster 
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NGC 1705-1). That is, we infer that the superbubble is not symmetrically located within the HI 
disk of NGC 1705, but must be on the near side. Thus, the blow-out is apparently one-sided at 
present (cf. MacLow, McCray, & Norman 1989). 

4.3. The Warm Photoionized Gas 

The column density of gas along our line-of-sight is dominated by the warm ionized gas. 
This implies that the shell of swept-up material (zone 5 in the superbubble) has been mostly 
photoionized by Lyman continuum radiation from the starburst. It is natural to identify the gas 
seen in absorption in the FUSE data with the emission-line nebula studied by Marlowe et al. 
(1995) and Meurer et al. (1992). 

The kinematics are consistent with this identification. Recall that the FUSE sight-line is 
heavily weighted towards material within about 3 arcsec of the bright super star cluster NGC 1705- 
1. The echelle spectra of this region analyzed by Marlowe et al. (1995) and Meurer et al. (1992) 
show double-peaked Ha emission-line profiles, with the blueshifted (redshifted) peak corresponding 
to emission from the front (back) side of the expanding superbubble. The measured LSR radial 
velocity of the blueshifted component is 560±20 km s~ x along the different position angles covered 
by their combined data sets. Within the uncertainties, this agrees with the mean velocity for the 
corresponding FUSE absorption-lines (569±10 km s _1 ). 

Since we now have a mass for the warm photoionized gas (§ 4.2), its Ha luminosity (Marlowe 
et al. 1995) can be used to directly determine that the mean electron density in this material is 
n e ~ 1 cm~ 3 . We emphasize that this makes no assumptions about the volume filling- factor of the 
gas. Taking T ~10 4 K (typical of photoionized gas), the corresponding thermal pressure is P/k 
~ 2xl0 4 K cm~ 3 . Can the starburst in NGC 1705 keep this gas photoionized? Taking n e ~ 1 
cm -3 , the measured radius of the superbubble (~ 500 pc), and the Lyman continuum luminosity 
of NGC 1705 (Q = 10 52 s _1 ), the Stromgren thickness of the photoionized supershell will be ~10 21 
cm and the column density of the ionized layer is then ~10 21 cm~ 2 . This is in good agreement 
with the column density of photoionized gas we infer from the FU SE data. 

We can compare the basic dynamical properties of the warm ionized gas in NGC 1705 to the 
predictions of the superbubble model. In convenient units, the radius and expansion speed of the 
superbubble are given by: 

r = 0.7Lli 5 echAQ n- 1/5 4 /5 kpc (1) 

« = 41i W M o< V V /5 (2) 
for an adiabatic superbubble inflated by a kinetic energy injection rate L me ch,40 (i n units of 
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10 erg s 1 ) into a uniform medium with nucleon density uq for a time tj (in units of 10 years). 

The age of the dominant super star cluster (NGC 1705-1) is well-constrained to be about 10 7 
years (de Mello, Leitherer, & Heckman 2000). This agrees with the dynamical age (t = 0.6r/v) of 
the superbubble (Marlowe et al. 1995; Meurer et al. 1992). The estimated kinetic energy injection 
rate due to supernovae is L mec /i,40 = 2 (Marlowe, Meurer, & Heckman 1999). Finally, the HI column 
in the center of NGC 1705 implies a mean nucleon density of no ~1 cm~ 3 for an ISM thickness of 
1 kpc. Our estimate of the swept-up mass of gas in the supershell (5xl0 7 M ) divided by volume 
enclosed by the supershell (~ 5 x 10 64 cm 3 ) likewise implies no ~1 cm~ 3 . Equations 1) and 2) 
above then predict a radius of 0.8 kpc and an expansion velocity of 50 km s _1 for the supershell. 
Given the overly-simplistic model (spherical symmetry, constant density) , we regard the agreement 
with the data as satisfactory. 



The neutral gas probed by FUSE is significantly less blueshifted than the ionized gas, and 
must therefore have a different origin. For this gas to remain neutral, it must be optically-thick to 
the Lyman continuum radiation from the starburst. This condition implies that a slab of gas with 
a column density logiV = 20.2 located inside the superbubble (and thus at a distance r < 500 pc 
from an ionizing source with a Lyman continuum luminosity of Q = 10 52 s _1 ) must have a density 
n > 8 cm~ 3 . Recall that the mean ISM density in NGC 1705 is roughly 1 cm -3 . Moreover, logiV 
= 20.2 and n > 8 cm -3 implies that the neutral absorbers must be very small (smaller than a few 
pc). 

Thus, only relatively small dense clouds would be neutral. On-average, the clouds are out- 
flowing, so most can not be in any undisturbed ambient medium in front of the superbubble's 
outer shell. We therefore identify the outflowing material with clouds in the superbubble's interior. 
These were presumably clouds in the ISM of NGC 1705 that were overtaken and engulfed as the 
superbubble propagated through the more tenuous inter-cloud medium. This would be a natural 
consequence of a multi-phase ISM (e.g., White & Long 1991). 

The dynamical model here is one of dense clouds exposed to the outflowing hot gas in zones 
2 and 3 of the superbubble. We would expect to see absorbing material that is injected from 
quiescent material at or near v sys , and which is then accelerated up to some terminal velocity as it 
is transported outward by the hot outflow. Following Heckman et al. (2000), an interstellar cloud 
with column density N, originally located a distance tq from a starburst, will be accelerated by a 
spherically-symmetric outflow that carries an outward momentum flux p up to a terminal velocity 
given by: 



4.4. 



The Neutral Gas 



vterm = 100 (p / 10 32 dynes) 1 ' 2 (r / '10 21 cm)^ 1/2 (N / 10 : 



,20 



cm 



2 )- 1/2 km/s 



(3) 
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We have chosen values for these parameters that are appropriate to the neutral gas in NGC 1705. 
The predicted velocity range of absorption relative to v S y S would be ~ to -100 km s" 1 . This agrees 
tolerably well with the velocities of the blueshifted neutral gas observed in our FUSE data. 

4.5. The Coronal Gas 

The detection of the interstellar OVI absorption-line (the first in a starburst galaxy) is the 
most important result in this paper. We will organize our discussion of the coronal gas in NGC 1705 
as follows. 

We begin by making some simple inferences about the basic physical properties of the gas 
that are independent of any specific model for the origin of the gas (§ 4.5.1). In particular, in 
§ 4.5.2 we point out that - independent of the detailed thermal/dynamical history of the gas - the 
observed OVI column density and line width imply that we are observing gas that has been heated 
to T > 3 x 10 5 K and which then cools radiatively. These model-independent results can then be 
used to constrain the radiative cooling rate from the coronal gas and show that it is small compared 
to the supernova heating rate. 

We then turn our attention to specific models for the OVI. We will not exhaustively consider all 
the possible alternative models (for example, conduction fronts associated with either the supershell 
fragments or the HI clouds discussed above). Instead we will describe the shortcomings of the 
standard superbubble model (§ 4.5.3), and then consider what we regard as an especially plausible 
alternative model in which the OVI arises in a hydrodynamical interaction between hot outrushing 
gas and the dense, cooler fragments of the superbubble shell during the "blow-out" phase. This 
idea is considered analytically in §4.5.4 in the context of models of turbulent mixing layers, and 
then using numerical hydrodynamical models in §4.5.5. 

4-5.1. Basic Physical Properties 

The combination of a detection of blueshifted OVIXW32 in absorption and an upper limit to 
corresponding redshifted emission allows us to place an upper bound on the density and pressure 
in the coronal gas (since we know the column density, and have an upper limit to the emission- 
measure). This calculation only assumes reasonable symmetry between the front and back sides of 
the outflow averaged over the kpc-scale projected FUSE aperture. 

The OVI ion is abundant only over a rather narrow temperature range (e.g. Sutherland 
& Dopita 1993), so we will assume T ~ 3xl0 5 K. The upper limit to the extinction-corrected 
OV 1X1032 flux in the FUSE aperture corresponds to an upper limit to the OVI intensity of 
-^1032 < 2.5xl0 4 photons cm~ 2 s _1 ster -1 . Following Shull & Slavin (1994), our observed OVI 
column density then leads to an upper limit to the electron density in the coronal gas of n e < 0.1 
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cm 3 . The corresponding upper limit to the pressure is P/k < 6xl0 4 K cm 3 . 

We have estimated a pressure in the warm photoionized gas of P/k ~ 2xl0 4 K cm" 3 . The 
pressure in the coronal-phase gas probed by the OVI ion should be the same, since zones 3 through 
5 in the adiabatic superbubble model are isobaric (MacLow & McCray 1988). We note that the 
isobaric condition could be invalidated by the presence of a dynamically-significant magnetic field (in 
which case the thermal pressure in the coronal gas should be higher than in the denser photoionized 
material). However, our model-independent upper bound of P/k < 6x 10 4 K cm -3 in the coronal gas 
is only a factor of three higher than the pressure implied by an assumption of isobaric conditions. For 
T ~ 3xl0 5 K, our estimated pressure then implies that the characteristic density in the coronal gas 
is n e ~3xl0~ 2 cm" 3 . For a metal abundance of 1/8 Solar in NGC 1705 (Heckman et al. 1998) the 
implied radiative cooling time from T ~ 3xl0 5 K is t coo \ = lxlO 6 years (interpolating between the 
collisional-ionization-equilibrium models for different metallicities in Sutherland & Dopita 1993). 
The assumption of collisional-ionization-equilibrium is probably reasonable for OVI, since the 
relevant recombination times at our assumed density and temperature (Nahar 1999) are roughly 
an order-of-magnitude less than the radiative cooling time. 

The estimated column density in the coronal gas in NGC 1705 is ~10 19 cm -2 , so the above 
density implies a characteristic thickness for the absorbing material of ~ 100 pc. If this path length 
is contributed by N total clouds, the sound crossing time of a cloud is ~ 1.2 x lO 6 ^" 1 years. Thus, 
the ratio of sound-crossing and radiative cooling times in a cloud is ~ 1.2 /N and the assumption 
that the coronal gas is in rough pressure-balance with its surroundings is therefore plausible. 

We have estimated that the total mass of coronal-phase gas is ~ 6 x 1O 5 M . The above cooling 
time then implies a cooling rate of M ~ 0.6 Mq per year. This can be compared to the average rate 
at which mass has been swept up by the superbubble over its lifetime: M ~ M/t ~ 5 x 1O 7 M /1O 7 
years ~ 5 M per year. The implied cooling luminosity is E = 3/2kTM / [imu ~ 10 39 erg s _1 . 
This is about 5% of the estimated rate at which the starburst supplies mechanical energy. 

Thus, radiative cooling associated with the coronal gas should not dominate the dynamical 
evolution of the outflow. We note that radiative cooling from the hotter gas detected in soft X-rays 
(Lx ~ 10 38 erg s _1 - Hensler et al. 1998) is significantly smaller than our estimate of E, and is 
therefore negligible. 

4-5.2. The Origin of the Observed Column Density 

Edgar & Chevalier (1986) have computed the expected column densities of various ions for the 
generic situation in which gas is heated to a temperature To ~ 10 6 K and then cools radiatively. 
The total column density of cooling gas is just given by N coo i = Nt coo i, where t coo \ is the radiative 
cooling time and TV is the rate of cooling per unit area. For a flow speed v, mass conservation in the 
cooling flow implies v = N/uq. Thus, the cooling column can also be written at N coo i = n§t coo iv. 
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It is important to note that N coo i is independent of density, since t coo \ oc rip 1 . Moreover, 
at coronal temperatures the cooling is dominated by metals, so N coo i oc t coo i oc Z^ 1 (where Z 
is the metallicity). Since Novi <x N coo iZ, the characteristic OVI column density in radiatively 
cooling coronal gas is essentially independent of density and metallicity and depends only on the 
value for v = N/n . Edgar & Chevalier (1986) calculate that N OV i tC ooi ~ 4 x 10 14 cm" 2 (v/100 
km/s). Calculations of OVI production in high-speed radiative shocks give similar values (Dopita 
& Sutherland 1996). This fiducial column density agrees well with our measured values for Novi 
and v in NGC 1705. 

This good agreement has two immediate implications. First, independent of the detailed 
dynamical and thermodynamical history of the OVI, this gas it is almost certainly the result of 
the radiative cooling of initially hotter gas. Second, we can now estimate the implied cooling rate. 

Using the Edgar & Chevalier (1986) models, our observed value logNovi = 14.3 implies that 
N/uq = 5 x 10 6 cm s -1 . In their models the gas cools from To = 10 6 K (and we note that 
significantly higher values for To in NGC 1705 are excluded by its low X-ray luminosity). Since we 
estimate P/k = 2xl0 4 K cm" 3 , it follows that n = 0.01 cm" 3 at T = 10 6 K and thus N = 5 x 10 4 
cm -2 s _1 . To calculate the implied rate at which gas is cooling (M = NAmu), we take a surface 
area A = 6 x 10 43 cm 2 for the NGC 1705 superbubble (see above), and find that M = 0.07M© per 
year. As the gas cools from 10 6 K its cooling luminosity is E = 3/2kTM / fimu ~ 10 39 erg s _1 . 
This is about 5% of the estimated rate at which the starburst supplies mechanical energy, which 
agrees well with the more naive estimate in § 4.5.1 above. 

The above estimates of the cooling rate depend upon our assumed pressure. However, the 
upper limit to the luminosity of OV/A1032 emission-line (§ 3.2.3) yields a pressure-independent 
upper limit to the cooling rates of M < 0.3 M Q per year and E < 4xl0 39 erg s _1 (Edgar & 
Chevalier 1986). These limits are consistent with the above estimates. 

4-5.3. The Failure of the Simple Superbubble Model 

In the simple superbubble model there are two plausible origins for the coronal gas. First, if 
the speed of outer shock driven into the ambient ISM is high enough, OVI ions will be abundant 
behind this shock (in zone 5). This does not appear to be feasible in NGC 1705, since the minimum 
required shock speed (post-shock temperature) is ~ 150 km s" 1 (3xl0 5 K) - see Shull & McKee 
(1979) and Dopita & Sutherland (1996). This is significantly higher than the expansion speed 
measured in for the shell in NGC 1705: ~ 53 km s _1 , corresponding to a post-shock temperature 
of only 40,000 K. 

Second, thermal conduction will transfer heat from zone 3 to zone 5 and create coronal-phase 
gas at the interface (zone 4). Weaver et al. (1977) predict an OVI column density in this material 
of 
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N OVI = 1 x 10- W 35 ^2mo4 /35 -- 2 (4) 

where Zq is the oxygen abundance relative to Solar. For Zq = 1/8, the predicted value for 
Novi is a factor of ~15 smaller than we measure. In this case, Novi is smaller than the value 
for radiatively cooled gas (§ 4.5.2) and has a direct dependence on the metallicity. This is because 
cooling in the conductively heated zone is dominated by adiabatic expansion losses (Weaver et al. 
1977). 

The relative velocity of OVI absorption in NGC 1705 is also inconsistent with the simple 
superbubble model, which predicts that the outflow velocity in this material (zone 4) will be 
substantially less than the outflow speed of the superbubble shell (zone 5) . Our data instead show 
that the OVI outflow speed is probably even larger than that of the shell (77±10 vs. 53±10 km 
s _1 - see above). 

We conclude that the simple superbubble model does not account for the observed properties 
of the OVI absorption-line in NGC 1705. 

4-5.4- Hydrodynamical Heating During Blow- Out 

We have argued in § 4.5.2 that the OVI arises in a flow of gas that has been heated to an 
initial temperature of To > 3 x 10 5 -fT and which is then cooling radiatively. As discussed in § 4.5.3, 
this heating/cooling process does not correspond to the outer shock in the superbubble (zone 5) 
because the observed expansion velocity of the superbubble is too slow to produce OVI behind 
this shock. In this section, we therefore consider a plausible alternative model. 

We have argued above that that the superbubble in NGC 1705 is in the process of breaking 
out of the ISM of NGC 1705. During this phase, the expansion speed of the superbubble shell 
accelerates and Rayleigh- Taylor instabilities will cause the shell (zone 5) to fragment. This allows 
the hot X-ray emitting gas in zone 4 to push its way out through the fragmented shell, and 
Kelvin-Helmholtz instabilities are believed to lead to the turbulent mixing of this hot outrushing 
gas with the cooler shell fragments (e.g. MacLow, McCray, & Norman 1989). Slavin, Shull, & 
Begelman (1993 hereafter SSB) have investigated the emission and absorption-lines produced by 
the intermediate-temperature (coronal phase) gas in these "turbulent mixing layers" (TMLs). 

We emphasize that the kinematics of the OVI absorbers in NGC 1705 are consistent with 
TMLs. The blueshift of the OVI relative to the shell material is expected as the hot gas rushes 
out through "cracks" in the shell. The magnitude of this blueshift follows from SSB: the velocity 
of the intermediate temperature gas in the TML relative to the cool gas from which it is created 
is given by v tmi ~ v hot {T coo i/T hot ) 1 / 2 , where v hot is the relative flow speed between the gas at T coo i 
and the outrushing hot gas at T^ ot . For the photoionized shell fragments that dominate the mass 
in NGC 1705 T coo \ ~ 10 4 K. Based on X-ray spectroscopy of dwarf starburst galaxies (della Ceca 
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et al. 1996,1997; Hensler et al. 1998), T hot = fewxlO 6 to 10 7 K. As this hot gas flows past the 
shell fragments, its maximum relative velocity will be roughly its sound speed (i.e. v^ot ~ 500 km 
s _1 ). Thus, for T coo i/Thot ~ 10 -3 , vt m i =16 km s -1 for Vhot = 500 km s _1 . This can be compared 
with NGC 1705 where the OVI absorption line is blueshifted by 24±10 km s _1 relative to the 
photoionized shell material. 

As pointed out by SSB, the ionic column density per TML is independent of the pressure 
in the gas, but is directly proportional to the relative velocity of the cool and hot gas (vh ot ). 
Also the OVI column will be appreciable only if the maximum temperature of the TML exceeds 
~ 3 x 10 5 K. Their models only extend up to Vhot = 100 km s" 1 and logT tm i = 5.5, and these predict 
OVI columns that are of-order 10 12 cm -2 per TML. These values could be adjusted upward in 
TMLs with larger values for Vhot an d logTtmi- In terms of the arguments in § 4.5.2, the low OVI 
columns in the TML models are partly attributable to the small relevant flow velocity in the TML: 
N /no = vtmi = Vhot(Tcooi /Thot) 1 ^ 2 ~ 10 km s -1 . It appears that of-order 10 2 TMLs per line-of-sight 
are needed in NGC 1705 (we see a 'sea" of shell fragments, each with its own TML). Proportionately 
fewer TMLs would be needed for larger flow velocities. 

4.6. Insights from Numerical Hydrodynamics 

As an initial attempt to investigate qualitatively the dynamics of gas in the complex transitional 
stage between superbubble and superwind, we have performed several 2-dimensional hydrodynam- 
ical simulations of a superbubble blowing out of a dwarf galaxy. We calculate the column densities 
and velocity structure of simulated lines of sight through these models, in particular the relative 
velocities of the warm and coronal phases that the FUSE observations of NGC 1705 show can not 
be explained by the classic superbubble model. 

The simulations were performed using the hydrodynamical code described in Strickland & 
Stevens (2000), with the initial conditions altered to roughly correspond to NGC 1705. These 
simulations were performed in cylindrical coordinates, assuming rotational symmetry around the 
z-axis (the minor axis of the galaxy). The hydrodynamical grid covers a region 2.5 kpc in radius 
by 2.5 kpc high along the z-axis with 500 by 500 equally sized cells. 

The initial ISM was set up in rotating hydrostatic equilibrium. We use a King model for the 
gravitational potential, producing a rotational velocity at r = 7 kpc of 65 km s — 1 . The ambient 
gas was given a temperature of 1.6 x 10 4 K, to approximate turbulent pressure support of the 
ISM, based on the ~ 15 km s _1 velocity dispersion of the HI gas. The resulting initial ISM density 
distribution has a peak number density (in the nucleus) of 2.5 cm~ 3 , a total column density through 
the galaxy along the minor axis of Ah = 3.5 x 10 21 cm~ 2 , a vertical scale height of 0.2 kpc, and a 
total mass of M gas = 8 x 10 7 M Q (within r = 2.5 kpc). These properties are all crudely similar to 
those observed in NGC 1705 (Meurer et al. 1998). 

We modeled the return of mass and mechanical energy into the ISM from starburst event as 
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a single instantaneous burst of star- formation, using the low metallicity models (Z = O.25Z ) of 
Leitherer & Heckman (1995). Appropriate amounts of mass and thermal energy were added at each 
computational time step to those cells within our assumed starburst region, a cylindrical region 
150 pc in radius and 60 pc thick vertically. 

Our fiducial model has a time-average mechanical luminosity of L m ech = 2 x 10 39 erg s _1 , 
somewhat weaker than the observed starburst NGC 1705, but chosen to give blow-out at t ~ 13 
Myr (approximately the observed dynamical age of the superbubble in NGC 1705). At blow-out 
the radius of superbubble is ~ 500 pc in the plane of the galaxy, and about twice as large along 
the minor axis. 

We assume an ISM metal abundance of 1/10 Solar for the purposes of calculating radiative 
cooling and column densities. We also performed other simulations to crudely explore parameter 
space: First, a simulation with a more powerful starburst (L mec ft = 7 x 10 39 erg s _1 ), which blows 
out earlier at t ~ 8 Myr but otherwise is very similar to the fiducial model in other properties. 
Second, a model using Solar abundance to investigate the scaling of column density with metal 
abundance. Third, a higher resolution model with twice the resolution in every dimension, to 
investigate the dependence of calculated properties on numerical resolution. 

An example of the model is shown in Figure 4, which plots the density and temperature of 
the gas in our fiducial low-metallicity model at t = 12.5 Myr. The OVI absorption seem in these 
simulations arises primarily in material at the interfaces between cool dense gas, e.g. in the shell of 
the superbubble or shell fragments after blow-out, and the hot (T > 10 6 K) gas filling the bubble 
volume (see Figure 4). In almost all cases these interfaces are numerically unresolved. 

We calculate the column density as a function of line-of-sight velocity profiles in each model 
at three different epochs (pre-, during, and post- blow-out, typically at t = 10, 13 and 18 Myr 
respectively), for four different lines of sight from the nuclear starburst (assuming galaxy inclination 
angles of i = 0, 30, 60 k, 80°), and for two different ionic species that probe different observed phases 
in superwinds (warm gas: Nil; coronal gas: OVI). Fractional ion abundances as a function of 
temperature were taken from Sutherland & Dopita (1993), which assume collisional ionization 
equilibrium. This a reasonable assumption in the case of OVI (see above). It is also a reasonable 
mock-up of Nil (in reality, photoionization by stars in NGC 1705 keeps the gas containing Nil 
near T ~ 10 4 , close to the temperature floor of 1.6 xlO 4 K in the simulations). 

Much of the UV light in NGC 1705 comes from the compact central star cluster, but a sig- 
nificant fraction comes from a more diffusely spread component (Meurer et al 1995). The width 
of 30" FUSE aperture used corresponds to a projected width of ~ 900 pc at the distance of NGC 
1705, so the observed spectrum is a sum over many slightly different lines of sight. To investigate 
the effect of this on the simulated absorption line profiles we explored two cases, which should 
bracket the real situation. In the first case we assumed all the light comes from a single point-like 
source at the center of the galaxy, and in the second case that the light comes from uniform surface 
brightness disk of radius 300 pc. We found very little difference in the simulated absorption line 



- 21 - 




200 400 600 800 200 400 600 800 200 400 600 800 

r (pc) r (pc) r (pc) 



Fig. 4. — Grey scale images of gas density and temperature in our fiducial hydrodynamical simu- 
lation of a dwarf starburst, shown as the superbubble blows out of the disk at t = 12.5 Myr. The 
three panels show (a) log number density (cm -3 ), (b) log gas temperature (K) and (c) the location 
of the OFI-absorbing material (gas with 5.2 < log T (K) < 5.8). 



(a) 



10' 



10 



12 



10 



13 







10 



14 



"OVI 



(cm"") 



10 



15 



10 



16 




-300 -200 -100 

<v 0V| > (km s _1 ) 



Fig. 5. — (a) Histograms of OVI column density, covering all lines of sight and all epochs in the 
three low metallicity numerical simulations (solid line) and the simulation using Solar metallicity 
(dashed line), (b) Corresponding histograms of the mean OVI absorption line velocity (weighted 
by the OVI column density). 
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profiles between the two methods, and so quote only results using the point-like light source. 

For any particular model, the column densities and kinematics of the gas vary strongly between 
different sight lines, due to the structural complexity that arises in 2-D simulations (see Suchkov et 
al. 1994, Strickland & Stevens 2000). The total column densities (in either OVI or Nil) do not 
change systematically with time in these simulations. We find that the mean and variance in total 
column density is the same before, during and after blow-out, considering all 12 simulated lines of 
sight at each specific epoch. 

As shown in Figure 5, we find OVI column densities in the range log Nov I ~ 13 ± 0.5 in the 
low metallicity models (Z = 0.1Z©), and log Novi ~ 14 ± 0.5 in the Solar metallicity model. Only 
about 10% of the lines-of-sight in the low metallicity models have OVI columns comparable to 
the observed value in NGC 1705. One reason for the failure of these simulations to reproduce the 
observed OVI column densities may be that there is insufficient numerical resolution to resolve the 
cooling zones between hot and cool gas, where the OVI-absorbing gas exists. Evidence to support 
this idea comes from comparing the predicted column densities in the fiducial model with those in 
the high resolution model. Although the predicted columns for low ionisation species always agree 
to within a factor 2, the predicted OVI column densities agree to within a factor 2 in only 50% of 
the sight lines, with the remaining sight lines differing by a factor of up to 30. The lack of thermal 
conduction in these models may also be partially responsible for the low OVI columns, although 
the arguments presented in § 4.5.3 suggest thermal conduction alone can not produce the observed 
OVI columns. Our simulations do show that the column density of lower ionization material is 
considerably higher than that of coronal material, and that this phase contains most of the mass 
in the outflow (consistent with our observations). 

The column- weighted mean velocities for the coronal gas along the different sight-lines generally 
lie in the range vovi ~ —100 to ~ km s -1 , although the full range covers vovi ~ —320 to ~ +50 
km s _1 (Figure 5). This gas is typically more blue shifted, and has a wider range of outflow 
velocities, than the cooler gas in our simulations. For example, the typical velocity range for Nil 
is vnh ~ —70 to ~ km s _1 and the full range seen only covers vnii ~ —110 to ~ km s _1 . The 
velocity range covered appears very similar in both low and Solar metallicity models, so it appears 
that the kinematics of the gas in these simulations is independent of the metal abundance. For 
OVI, the mean outflow velocity (and range of velocities), does appear to increase with time, moving 
from pre-blow-out through to post-blow-out. This is not the case for lower ionization species, e.g. 
the distribution of mean velocity for Nil is very similar at all three epochs. 

In summary, while the numerical models underpredict the OVI columns, they can explain 
one aspect of the observed data on NGC 1705. The mean outflow velocity of high-ionization gas 
probed by OVI is as high as, or higher than, that of low-ionization gas. This is in contrast with 
the standard superbubble model in which the outflow speed in the coronal gas (zone 4) is less than 
that of the superbubble shell. This result supports the idea that blow-out is responsible for the 
kinematics of the coronal phase in NGC 1705. 
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5. Implications 

As summarized in § 1, probably the greatest uncertainty in assessing the impact of starburst- 
driven outflows on the evolution of dwarf galaxies is the importance of radiative cooling (which 
will directly determine the outflow's dynamical evolution). In the absence of severe radiative 
cooling, both analytic arguments and hydrodynamical simulations imply that outflows like the one 
in NGC 1705 will easily blow-out of their galaxy's ISM (e.g., MacLow & McCray 1988; Marlowe 
et al. 1995; Martin 1998,1999; MacLow & Ferrara 1999; Strickland & Stevens 2000). The hot gas 
that blows out of the superbubble will then be able to heat and chemically enrich the inter-galactic 
medium (e.g., Ponman et al. 1999; Tozzi, Scharf, Sz Norman 2000) and/or an extended gaseous 
halo surrounding the dwarf galaxy (Silich & Tenorio-Tagle 1998). 

We will not repeat these arguments here, but simply note that within the context of the 
standard McKee & Ostriker (1977) formalism, the lack of radiative cooling would be a result of 
the exceedingly high supernova rate per unit volume, which leads to a porosity near unity in the 
ISM (i.e., the ISM within which the supernovae detonate is filled almost entirely with hot, diffuse 
gas). This is especially pertinent to NGC 1705, in which nearly half of the supernovae will detonate 
inside the pc-scale super-star-cluster NGC 1705-1. 

Existing X-ray data establish that radiative cooling from hot (T > 10 6 K) gas is energetically 
insignificant in NGC 1705 (Hensler et al. 1998) and other starburst galaxies (cf. Heckman 2000 
and references therein). The new data from FUSE are very instructive, since they directly probe 
for the coronal-phase (T ~ 10 5 - 10 6 K) gas near the peak of the cooling curve where hitherto- 
unobserved radiative cooling might be severe. As we showed in § 4.5.2 above, the observed OVI 
column density (1.8 x 10 14 cm~ 2 ) and flow speed (~ 10 2 km s" 1 ) implies that we are seeing gas in 
the outflow that has been heated to T > 3 x 10 5 K and then radiatively cooled. We have estimated 
that the implied rate of radiative cooling is only ~6% of the rate of supernova heating (§ 4.5.2). 
Thus, there is no place left to hide the radiation necessary to quench the supernova-driven outflow 
in NGC 1705. 

What then is the fate of the outflowing gas? For an isothermal gravitational potential that 
extends to a maximum radius r max , and has a virial velocity v ro t, the escape velocity at a radius r 
is given by: 

Vesc = v rot [2(l + ln(r max /r))) l/2 (5) 

The rotation curve for NGC 1705 implies v ro t = 62 km s _1 (Meurer et al. 1998), and the 
rotation curve remains relatively flat out to the last measured point at r = 7 kpc. Thus, the 
minimum escape velocity at the present location of the superbubble wall (r ~ 0.5 kpc) will be 
roughly 170 km s -1 . This is significantly higher than the velocities of most of the outflowing gas 
that we observe in-absorption in NGC 1705. Thus, there is no direct evidence that most of this 
material will escape NGC 1705. Maps of the HI in NGC 1705 (Meurer et al. 1998) show that the 
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large-scale ISM in NGC 1705 will not be strongly affected by the outflow (which will only puncture 
a few-kpc-scale hole in the central region of the HI disk). 

On the other hand, the bulk of kinetic energy and most of the newly-created metals will reside 
in the hotter gas that is driving the superbubble's expansion. The temperature of this gas will be 
> 10 6 K (della Ceca et al. 1996,1997; Hensler et al. 1998) . We have argued that we are witnessing 
the blow-out phase in NGC 1705 in which this gas is vented through the fragmented superbubble 
wall. Following blow-out, in the absence of significant radiative cooling this hot gas will feed a 
galactic wind with a terminal velocity of at least 300 km s _1 (Chevalier & Clegg 1985). This is 
comfortably above v esc for any plausible value for r max . Once the hot gas vents out of the ruptured 
superbubble, there is no known external medium to impede the resulting high-speed wind. 

We conclude that NGC 1705 can eject most of the newly-created metals and kinetic energy 
returned by the central starburst over the past ~ 10 7 years. If typical of dwarf starbursts, this 
process could help create the low metal abundances in dwarf galaxies (e.g., Dekel &; Silk 1986; 
Lyden-Bell 1992) and contribute to the heating and chemical enrichment of the intergalactic medium 
(Ponman et al. 1999; Gibson et al. 1997; Ellison et al. 2000). 

6. Conclusions 

We have presented new FUSE far-UV spectroscopy of the prototypical dwarf starburst galaxy 
NGC 1705. Previous optical and UV spectroscopy has established that the starburst in this galaxy 
is driving a large-scale outflow (Meurer et al. 1992; Marlowe et al. 1995). The FUSE data are es- 
pecially important because they probe the coronal-phase (T ~ 10 5 to 10 6 K) gas that may dominate 
the radiative cooling of the outflow, and thereby largely determine its dynamical evolution. The 
high quality of the data (spectral resolution of ~ 30 km s _1 and a S/N ~ 16) gives us important 
new insight into all the phases in the ISM of this galaxy. 

Firstly, we do not detect any H2, with an upper limit of logN}j2 < 14.6. This is similar to 
the results obtained by Vidal-Madjar et al. (2000) for the dwarf starburst IZwl8. In the specific 
case of NGC 1705 the low H2 abundance is consistent with Galactic sight-lines with similarly low 
E(B - V) (Savage et al. 1977). 

We do detect absorption from three other phases of the ISM: neutral gas, warm photoionized 
gas, and coronal-phase gas (probed with the OKTA1032 line). The total column densities (cm -2 ) 
in the three phases are logN = 20.2, 20.9, and >19.0 respectively. The first is directly measured 
from Lyman series, while the latter two are estimated from the measured ionic columns assuming 
a metal abundance of 12% solar (as measured in the optical emission-line gas). 

All the interstellar absorption-lines associated with NGC 1705 are broad (~ 10 2 km s _1 ) and 
blueshifted. Thus, gas is flowing out of NGC 1705 with (v — v sys ) = -32, -53, and -77 km s _1 for 
the neutral, photoionized, and coronal gas respectively. 
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The mass and kinetic energy in the outflowing gas probed by FUSE is dominated by the warm 
photoionized gas, which is also seen via its optical line-emission. The size, morphology, outflow 
speed, mass, kinetic energy, and dynamical age of this structure are all consistent with a simple 
model of an adiabatic superbubble whose expansion is driven by a piston of hot gas created by the 
cumulative effect of the supernovae in the starburst (e.g., Weaver et al. 1977). 

In contrast, we have shown that the neutral gas must reside in relatively small dense clumps 
in order to be optically-thick in the Lyman continuum. We attribute the outflowing material to 
clouds in the ISM that have been overtaken by the superbubble and accelerated by the outflowing 
hot gas. 

The properties of the coronal gas probed with 0VI\1Q2>2 are inconsistent with the predictions 
of the simple superbubble model. The observed expansion speed of the superbubble (53±10 km 
s _1 ) is too small to produce OVI behind its shock-front. The blueshift of the OVI line is larger 
than that of superbubble shell (while the model predict that it should be substantially smaller). 
Finally, the observed OVI column (logNovi = 14.3) is a factor of ~15 too large compared to 
the superbubble model (in which such gas arises in a conductively-heated and adiabatically-cooled 
interface between hot gas and the cool outer superbubble shell). 

We therefore proposed the following origin for the OVI absorber, and have investigated its 
plausibility via numerical hydrodynamical models. The morphology and relative size-scale of the 
superbubble implies that it has begun to "break out" or "blow out" of the ISM in the disk of 
this galaxy. During the blow-out phase, the superbubble shell will accelerate and then fragment 
(e.g., MacLow & McCray 1988). As it fragments, the piston of hot gas inside the superbubble 
will flow out through the "cracks", and the resulting hydrodynamical interaction between this 
outrushing gas and the shell fragments will create intermediate-temperature coronal gas that can 
produce the observed OVI absorption (e.g., Slavin et al. 1993). This process accounts nicely for 
the observed kinematics of the OVI absorption. We also emphasized that for the observed flow 
speed of ~ 10 2 km s _1 , the observed OVI column density is just what is expected for gas that has 
been collisionally-ionized and then radiatively-cooled, independent of its density or metallicity and 
independent of any specific hydrodynamical model (Edgar & Chevalier 1986). 

We have argued that the coronal gas should be in rough pressure-balance with the warm 
photoionized gas (whose pressure we determined to be P/k = 2 x 10 4 K cm -3 ). This allowed us 
to estimate cooling rates of ~0.07 M & per year and ~10 39 erg s _1 in the coronal gas, based on 
the OVI absorption-line. Independent of the gas pressure, the lack of redshifted OVI emission 
from the backside of the outflow implies upper limits to the cooling rates of <0.3 M Q per year and 
< 4 x 10 39 erg s _1 . The rate of radiative cooling due to hotter (T > 10 6 K) gas is negligible in 
comparison (Hensler et al. 1998). Since the rate of supernova-heating in NGC 1705 is 2 x 10 40 erg 
s _1 (Marlowe et al. 1999), we have argued that radiative cooling is not dynamically important in 
the NGC 1705 outflow. 

In the absence of significant radiative cooling, the superbubble in NGC 1705 should be able to 
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blow out of the ISM. Based on the relatively low outflow velocities, it is unlikely that the absorption- 
line gas we observe with FUSE will escape from the galaxy. Likewise, the bulk of the global ISM 
(HI) in NGC 1705 will be retained (cf. De Young k Heckman 1994; MacLow k Ferrara 1999). 
However, we have argued that the hotter gas that drives the outflow (and which is now being 
vented out of the ruptured superbubble) can escape and thereby carry away most of the kinetic 
energy and newly-created metals supplied by the starburst. This process has potentially important 
implications for the evolution of both dwarf galaxies and the IGM. 

We thank the members of the FUSE team for providing this superb facility to the astronomical 
community. We thank Colin Norman and Robin Shelton for enlightening discussions, and an 
anonymous referee for a detailed and thoughtful report that improved the paper. This work was 
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through Chandra Postdoctoral Fellowship Award Number PFO-10012, issued by the Chandra X-ray 
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Table 1. Equivalent Widths and Radial Velocities 



ID 


A 


Wmw 


WjV1705 


VMW 


V JV1705 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


H2 Q(l) 


1009.770 


<35(3<t) 


<35(3<t) 






C II 


1036.337 


510+30 (LI) 
492±22 (L2) 


>680 (Ll)c 
>670 (L2)c 


+18 


+582 


C III 


977.020 


560±30 (S2) 
520+35 (bl) 


>800 (S2)b 

^_ rrriA /n-| \ l 

>700 (bl)b 


+28 


+555 


N I 


1134.165 


121±10 (LI) 

i n a 1 1 n /T ri\ 

134+10 (L2) 


23+10 (LI) 
24+10 (L2) 


+28 


+582 


N I 


1134.415 


155±10 (LI) 
156±10 (L2) 


32+10 (LI) 
23+10 (L2) 


+24 


+590 


N I 


1134.980 


171+11 (LI) 
161±10 (L2) 


35±11 (LI) 
35±10 (L2) 


+25 


+602 


N II 


1083.990 


356+25 (S2) 


338±35 (S2) 


+7 


+573 


O I 


1039.230 


260+20 (Ll)a 
252+20 (L2)a 


243±15 (LI) 
237±15 (L2) 


+14 


+592 


O I 


988.773 


302+25 (b2)a 
283±30 (SI) 


480±30 (S2) 
540±45 (SI) 


+10 


+595 


O VI 


1031.926 


252±22 (LI) 
220±22 (L2) 


171±18 (LI) 
191±19 (L2) 


+25 


+545 


Si II 


1020.699 


129±16 (LI) 
104+12 (L2) 


94±24 (LI) 
65±20 (L2) 


+ 16 


+589 


S III 


1012.502 


143+20 (Ll)b 
128±20 (L2)b 


250±15 (LI) 
215±20 (L2) 


+8 


+578 


S IV 


1062.662 


<36 (3cr) 
<36 (3a-) 


85±15 (LI) 
82±15 (L2) 




+555 


Ar I 


1048.220 


76+10 (LI) 
83+10 (L2) 


60±10 (LI) 
80±20 (L2) 


+20 


+566 


Fe II 


1063.176 


125±11 (LI) 
140±13 (L2) 


162±16 (LI) 
177±21 (L2) 


+11 


+597 


Fe II 


1144.938 


192±12 (LI) 
206±11 (L2) 


262±17 (LI) 
241±17 (L2) 


+20 


+601 



Note. — The line equivalent widths in the Milky Way (Col. 3) and 
NGC 1705 (Col. 4) are in mA. These widths and the la errors were de- 
rived according to the prescription outlined by Sembach & Savage (1992). 
The values account for statistical noise and modest continuum placement 
uncertainties. In some cases, additional uncertainties (not listed) due to 
stellar blending and fixed-pattern noise introduced by the FUSE detectors 
may be warranted. The detector segment utilized (LiFl, LiF2, SiCl, and 
SiC2) is indicated in parentheses following each equivalent width. The let- 
ters appended after these measurements have the following meaning: a - 
Value may be affected by terrestrial O I airglow emission; b - Line probably 
has some stellar blending, measurement uncertain; c - Lower limit because 
of blending with other lines. The radial velocities of the centroids for Milky 
Way (Col. 5) and NGC 1705 lines (Col. 6) are in the LSR frame. These were 
derived by fitting single Gaussian profiles to the absorption lines. Typical 
uncertainties are ±10 km s . The wavelength scale in our data has been 
adjusted so that the Milky Way neutrals agree with the mean velocity of the 
lines from similar species in our HST STIS echelle spectrum of NGC 1705 
(Sembach et al. 2001; Heckman et al. 2001). 
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Table 2. NGC 1705 Column Densities 



Species 


A 


log(fX) 


logN 


logN 


logN 








(Sidel) 


(Side2) 


Final 


H2 (J=0) 






<14.2 


<14.2 


<14.2 


H2 (J=l) 






<14.2 


<14.2 


<14.2 


H2 (J=2) 






<14.2 


<14.2 


<14.2 


H I 


1025.7 


1.909 


20.2 


20.2 


20.2±0.2 


N I 


1134.1 


1.182 


<14.17 


<14.30 


<14.30 




1134.4 


1.483 


<14.08 


<14.07 


<14.07 




1135.0 


1.660 


13.97 


13.98 


13.97±0.08 


N II 


1084.0 


2.048 




14.70 


14.70±0.15 


I 


1039.2 


0.980 


15.64 


15.63 


15.63±0.08 


VI 


1031.9 


2.137 


14.26 


14.26 


14.26±0.08 


Si II 


1020.7 


1.460 


14.68 


14.51 


14.61±0.14 


S III 


1012.4 


1.556 


15.02 


15.00 


15.02±0.12 


S IV 


1062.7 


1.628 


14.45 


14.43 


14.45±0.08 


Ar I 


1048.2 


2.408 


13.44 


13.49 


13.46±0.10 


Fe II 


1144.9 


2.084 


14.53 


14.54 


14.54±0.10 




1063.2 


1.805 


14.54 


14.55 


14.54±0.10 



Note. - The /-values are taken from Morton (1991), and 
A is in A. Column densities (cm -2 ) were derived by converting 
the observed absorption profiles into optical depth profiles and 
integrating over velocity (see Savage & Sembach 1991). These 
column densities do not have any saturation corrections applied, 
though these are expected to be small given the great breadth 
of the lines. A comparison of the results for the two Fell lines 
listed indicates that this is indeed the case. In the case of the 
HI column, we quote the average of the value obtained from 
the higher-order Lyman series in the FU SE data and the value 
derived from Lya by Heckman & Leitherer (1997). 



